A Measurement of Atmospheric Neutrino Flux Consistent with Tau Neutrino Appearance 
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A search for the appearance of tau neutrinos from z/ M <-> v T oscillations in the atmospheric neutrinos has 
been performed using 1489.2 days of atmospheric neutrino data from the Super-Kamiokande-I experiment. A 
best fit tau neutrino appearance signal of 138 ± 48 (stat.) ig 2 (sys.) events is obtained with an expectation of 
78 ± 26 (sys.). The hypothesis of no tau neutrino appearance is disfavored by 2.4 sigma. 
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Atmospheric neutrino oscillations have been observed by 
several experiments QJ 0, S HI ; m particular, Super- 
Kamiokande (Super-K) has reported the first evidence for the 
sinusoidal signature of muon neutrino disappearanceJ16j] and 
made measurements of the oscillation parameters |7]. The 
Super-K atmospheric neutrino data favor <-> v T oscilla- 
tions and have excluded <-> v e ||l|] and pure v^, «-> ^sterile 
oscillations as a dominant source of the deficit of muon 
neutrinos. All of these results rely largely on the disappear- 
ance of muon neutrinos from the atmospheric neutrino flux, 
with no explicit observation of the appearance of tau neutri- 
nos and their charged-current (CC) weak interactions. In this 
letter, we analyze the Super-K atmospheric neutrino data in 
an attempt to demonstrate the appearance of v T interactions in 
the detector. 

Detection of CC v T interactions in the atmospheric neu- 
trinos is challenging for two reasons. First, the neutrino en- 
ergy threshold for tau lepton production is 3.5 GeV. The at- 
mospheric neutrino flux above this energy is relatively low. 
Assuming two flavor maximal mixing of «-> v T with Am 2 
= 2.4 x 10 -3 eV 2 , approximately one CC v T event is expected 
to occur in an atmospheric neutrino detector per kiloton-year 
of exposure. This corresponds to an estimated total of 78 v T 
events in the data sample presented. Second, the tau lepton 
has a short lifetime (290 femtoseconds) and decays immedi- 
ately into many different final states. These final states consist 
of electrons, muons, or one or more pions (plus always a tau 
neutrino). The recoiling hadronic system may also produce 
multiple particles. Water Cherenkov detectors such as Super- 
K are not suited for identifying individual CC v T interactions 
as there are generally multiple Cherenkov rings with no eas- 
ily identified leading lepton. Thus, we employ likelihood and 
neural network techniques to discriminate tau neutrino events 
from atmospheric neutrino events on a statistical basis. 

Super-Kamiokande is a 50-kton water Cherenkov detector, 
with a rock overburden of 2700 m water equivalent, located in 
Kamioka, Gifu prefecture in Japan. The detector consists of 
two concentric optically separated detector regions; the inner 
detector (ID) instrumented with 1 1,146 inward facing 20 inch 
diameter photomultiplier tubes (PMT) and the outer detector 
(OD) instrumented with 1,885 outward facing 8 inch PMTs. 
The details of the detector, calibrations, data reduction, and 
detector simulation can be found in Refs. flU. 

In this letter, the atmospheric neutrino data accumulated 
during the Super-K-I period (from May, 1996 to July, 2001; 
1489.2 live-day exposure) are analyzed. The atmospheric neu- 
trino events in Super-K are classified as fully-contained (FC), 
partially-contained (PC), and upward-going muon events |7]. 
In the present analysis, only FC events are used. FC events 
deposit all of their Cherenkov light inside the ID, from which 
the direction and the momentum of charged particles are re- 
constructed. The particle type is identified as "e-like (show- 
ering)" or "/i-like (non-showering)" for each Cherenkov ring 
based on the light-pattern |7|. 

Both tau neutrino and atmospheric neutrino (y e and v^) in- 
teractions in Super-K are modeled using a Monte Carlo (MC) 



simulation, NEUT ItHigIi. In our v T MC, only CC v r inter- 
actions are simulated. CC v T interactions are mostly deep- 
inelastic scattering (approximately 63%) due to the high en- 
ergy threshold for tau lepton production. The decays of tau 
leptons are simulated by the tau decay library, TAUOLA (Ver- 
sion 2.6) 1 12]. The polarization of tau leptons produced via 
CC v T interactions is also implemented from calculations by 
Hagiwara et al. 

The event signatures of tau neutrinos are characterized by 
the decays of tau leptons produced in CC v T interactions. In 
this analysis, we concentrate on the hadronic decays of tau 
leptons (approximately 65% of tau lepton decays). The shape 
of events containing the hadronic decays of tau leptons has a 
more spherical topology than that of backgrounds. Also, the 
extra pions produced in tau lepton decays can be tagged by 
looking for their decays and ring signatures in the Super-K 
detector. The primary backgrounds for v T signals are atmo- 
spheric neutrinos producing multiple pions via deep-inelastic 
scattering interactions. The following v T event selection cri- 
teria are applied to reduce the backgrounds: (1) The vertex 
must be reconstructed in the fiducial volume (2 m from the 
ID PMT surface) with no activity in the OD region, (2) Visi- 
ble energy {E V i S ) must be greater than 1.33 GeV, and (3) The 
most energetic ring must be e-like. These criteria can reduce 
the backgrounds by approximately 90% since the neutrino en- 
ergy threshold in CC v T interactions is higher than in most of 
atmospheric neutrino interactions and tau lepton decays have 
a high average multiplicity (resulting primarily in hadronic 
shower, i.e. e-like events). The effects of the event selection 
criteria are summarized in the top 3 lines of Table U 

In accordance with the event shape and characteristics of 
tau lepton decays, we define a set of five variables to further 
discriminate v T signals from backgrounds, which are: (a) Vis- 
ible energy, (b) Maximum distance between the primary inter- 
action and electron vertices from pion and then muon decays, 
(c) Number of ring candidates, (d) Sphericity in the laboratory 
frame, and (e) Clustered sphericity in the center of mass frame 
(FIG. DJ. The first two variables are Super-K standard vari- 
ables 1 7]. The number of ring candidates is determined with 
a ring-finding algorithm, which is sensitive to ring-fragments, 
and the last two variables, sphericity and clustered sphericity, 
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Generated in fiducial volume 




17135(100%) 


78.4 (100%) 


E ms > 1.33 GeV 


2888 


2943 (17.2%) 


51.5 (65.7%) 


Most energetic ring e-like 


1803 


1765 (10.3%) 


47.1 (60.1%) 


Likelihood > 


649 


647 (3.8%) 


33.8 (43.1%) 


NN output > 0.5 


603 


577 (3.4%) 


30.6 (39.0%) 



TABLE I: Summary for the numbers of events in data, atmospheric 
neutrino (f e ,/i) background MC, and tau neutrino MC after apply- 
ing each of v T event selection criteria and either the likelihood or 
the neural network cut, which is applied separately for each analysis. 
Neutrino oscillation is considered in the MCs with the oscillation pa- 
rameters: Am 2 = 2.4 x 10 -3 eV 2 and sin 2 28 = 1 .0, and the numbers 
are normalized by the live time of the data sample. 



3 



150 
100 
50 


w 

c 80 

<D 

W 60 
o 

S 40 

-O 

£ 20 
Z 





(a) 


■ . , 1 , 




3.5 


log(E vis ) log([MeV]) 


- 1 1 


i, (0 











10 20 
Ring Candidates 



30 




-1.5 -1 -0.5 
log(Sphericity) 




0.25 0.5 0.75 1 
Clustered Sphericity 



FIG. 1: The distributions of variables: (a) Visible energy, (b) Maxi- 
mum distance between the primary interaction and electron vertices 
from pions and then n decays, (c) Number of ring candidates, (d) 
Sphericity in the laboratory frame, and (e) Clustered sphericity in 
the center of mass, after applying the v T event selection criteria for 
downward-going data (points), v T MC (shaded histogram), and at- 
mospheric i/ e , M background MC events (solid histogram). (The his- 
tograms of v T MC are normalized arbitrarily.) In the likelihood anal- 
ysis, the data sample is divided into 5 energy bins. 
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FIG. 2: The likelihood (top) and NN output (bottom) distributions 
of downward-going (left) and upward-going (right) events for data 
(points), atmospheric neutrino background MC (dashed histogram), 
and the best fit including tau neutrino and backgrounds (solid his- 
togram). The shaded area shows a fitted excess of tau neutrino events 
in the upward-going direction. The events for likelihood C > or 
NN output > 0.5 are defined to be tau-like. 



are derived from a jet-based energy flow analysis. 

The energy flow analysis uses Cherenkov patterns to de- 
convolve the measured light distribution in Super-K. By asso- 
ciating the deconvoluted power spectrum with pseudo parti- 
cles, "jets" are reconstructed, and event shape variables such 
as sphericity are obtained. Sphericity (0 < S < 1) measures 
the spherical symmetry of an event and has a quadratic mo- 
mentum dependence giving more weight to higher momentum 
particles Il4ll5ll . Also, it is not invariant under the Lorentz 
transformation. Therefore, the sphericity calculated with dif- 
ferent clustering and reference frames probes different event 
characteristics. 

The expected distributions of each variable after apply- 
ing the v T event selection criteria for both signals and back- 
grounds are plotted in FIG. ^ The background MC is com- 
pared with downward-going data, where no v T appearance 
signal is expected as the probability for oscillating into v T 
is very small for the given path length and the measured at- 
mospheric Am 2 . The agreement of the downward-going data 
and background MC indicates the variables chosen for this 
analysis are well modeled by our MC simulation. The small 
discrepancy in log(Sphericity) between the data and the atmo- 
spheric neutrino MC is consistent with the systematic uncer- 
tainty in the vertex position. 

We have constructed a likelihood function using the five 
variables described above. The data sample is divided into 5 



energy bins: (1) E ms < 2.0, (2) 2.0 < E ms < 3.0, (3) 3.0 
< E ms < 6.0, (4) 6.0 < E ms < 12.0, (5) 12.0 < E ms [GeV]. 
The likelihood distributions for downward-going and upward- 
going events are shown in FIG. |2] The events for likelihood 
L > are defined to be tau-like. The likelihood distributions 
of data and background MC events agree for downward-going 
events. The agreement validates our analysis method. TableU 
summarizes the number of data, atmospheric neutrino back- 
ground MC, and tau neutrino MC events after applying all of 
the event selection criteria. 

A neural network (NN) is also trained with the five vari- 
ables. The network has 6 input neurons, 10 hidden neu- 
rons, and one sigmoid output neuron and is trained using 
back-gropagationby use of the MLPFIT neural network pack- 
age m. 

The distributions of the NN outputs are shown in FIG. |2] 
NN output > 0.5 is defined to be tau-like. The discrepancy 
between the data and MC just below the NN cut (NN output = 
0.5) is consistent with the systematic uncertainty in the deep- 
inelastic scattering cross section near the threshold, which is 
included in the overall systematic uncertainty estimation de- 
scribed below. 

After selecting the tau-enriched sample by applying the v T 
event selection criteria with either the likelihood (C > 0) or 
the neural network (NN output > 0.5) cut, the zenith angle dis- 
tribution is fitted with a combination of the expected tau neu- 
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FIG. 3: The zenith angle distributions for the likelihood (top) and 
neural network (bottom) analyses. Zenith angle cos# = — 1 (cos# 
= +1) indicates upward-going (downward-going). The data sample 
is fitted after v r event selection criteria are applied. The solid his- 
togram shows the best fit including v T , and the dashed histogram 
shows the backgrounds from atmospheric neutrinos (v a and v^). 
A fitted excess of tau-like events in the upward-going direction is 
shown in the shaded area. 



trino signals resulting from oscillations and the predicted at- 
mospheric neutrino background events including oscillations. 
The fitted zenith angle distribution and the x 2 , which is mini- 
mized, are 



N t otal(cOS0) = aiVtau + /3iV bkg , 
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where 7V° bs is the number of the observed events, A^ tau is 
the number of predicted tau neutrino events, A^ bkg is the MC 
predicted number of atmospheric neutrino background events, 
and er, is the statistical error for the i-th bin. The sample nor- 
malizations, a and (3, are allowed to vary freely, and the zenith 
angle distribution is divided into 10 bins, from —1 to 1 (cos# 
= — 1 (cos6> = 1) refers to upward-going (downward-going) 
events). The results of the fit are shown in FIG.|3] 

The minimum x 2 f° r the likelihood fit (the NN fit) is 
xL« = 7 - 6 / 8 D0F (9-8/8 DOF), and the X 2 assuming no 
tau neutrino appearance is 16.3/9 DOF (18.2/9 DOF). The 
normalizations of the best fit for the likelihood fit (the NN fit) 
are a = 1.76 (1.71) and (i = 0.90 (0.99). After correcting for 
efficiencies, these correspond to a best fit tau neutrino appear- 
ance signal of 138 ±48 (stat.) (134±48 (stat.)) for the like- 
lihood analysis (the NN analysis). As can be seen in FIG.|3] 



Systematic uncertainties for expected v T 


LH (%) 


NN (%) 


Super-K atmospheric v oscillation analysis 


21.6 


20.2 


(23 error terms) 






Tau related: 






Tau neutrino cross section 


25.0 


25.0 


Tau lepton polarization 


7.2 


11.8 


Tau neutrino selection efficiency 


0.4 


0.5 


LH selection efficiency 


4.8 




NN selection efficiency 




3.0 


Total: 


32.6 


34.4 




Systematic uncertainties for observed v T 


LH (%) 


NN (%) 


Super-K atmospheric v oscillation analysis: 






Flux up/down ratio 


6.5 


5.7 


Flux horizontal/vertical ratio 


3.6 


3.2 


Flux K/-7T ratio 


2.4 


2.8 


NC/CC ratio 


4.3 


3.8 


Up/down asym. from energy calib. 


1.4 


< 0.1 


Oscillation parameters: 






0.0020 < Amis < 0.0027 eV 2 


+5.8 


+8.8 




-2.6 


-3.3 


0.93 < sin 2 2023 < 1.00 


-3.3 


-3.9 


0.0 < sin 2 20i3 < 0.15 


-20.6 


-17.9 


Total: 


+ 10.7 


+ 12.0 




-22.9 


-20.3 



TABLE II: Summary of systematic uncertainties for the expected 
number of v T events (top) and for the observed number of v T events 
(bottom). The best fit values for each error term are listed for both 
likelihood (LH) and neural network (NN) analyses. 



an excess of tau neutrino signals is observed in the upward- 
going direction, and the data distribution agrees better with 
the prediction including tau neutrino appearance estimated by 
MC. The backgrounds that remain after applying all of the 
v T event selection criteria are mostly deep-inelastic scattering 
(CC DIS: 61.4% and NC DIS: 27.1%). 

Approximately 82.9 + 3.0% of events are in common to 
the tau-enriched samples selected by both analyses, for which 
Monte Carlo predicts 83.1% of events overlap. The results for 
the likelihood and the neural network analyses are consistent. 

The systematic errors from Super-K atmospheric neutrino 
oscillation analysis are re-evaluated for the present analysis; 
however in this estimation, the uncertainty in the absolute nor- 
malization is assumed to be 20%. All error terms except for 
those affecting Sub-GeV, PC, and upward-going muon events 
are considered. A detailed description of these uncertainties 
can be found in Ref . 1 7 ] . In addition, the uncertainties related 
only to the present analysis such as the v T cross section, v T 
polarization, tau likelihood, etc. are considered. The system- 
atic uncertainties for the expected number of v T events are 
summarized in TablellTI 

In determining the systematic uncertainties for the observed 
number of v r events, various effects (such as up/down ra- 
tio) that could change the up-down asymmetry of the back- 
ground MC and the data are considered. The systematic er- 
rors due to uncertainties in the oscillation parameters, Am| 3 
and sin 2023, are also estimated by using 68% C.L. allowed 
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parameter region obtained by the L/E analysis from Super- 
K 1 6]. The uncertainty due to knowledge of the sin 2 29 13 
is estimated using the limit obtained by the CHOOZ reactor 
neutrino experiment [17]. This systematic error is asymmet- 
ric because Multi-GeV electrons are expected to appear in the 
upward-going directions for non-zero #13, which would be 
backgrounds for v T signals. Table ITU shows the summary of 
systematic uncertainties. We also performed a study to check 
dependency on Monte Carlo neutrino interaction models using 
another model, NUANCE 11811 . The difference in the results 
is negligible. 

Combining these errors with the fit result, we obtain a best 
fit tau neutrino appearance signal of 138 ±48 (stat.) (sys.) 
from the likelihood analysis, which disfavors the no tau neu- 
trino appearance hypothesis by 2.4 sigma. This is con- 
sistent with the expected number of tau neutrino events, 
78 ±26 (sys.) for Am 2 = 2.4 x lCT 3 eV 2 , assuming the full 
mixing in «-> v r oscillations. 

In conclusion, the search for the appearance of tau neutri- 
nos from i/f_i <-> v T oscillations in the atmospheric neutrinos 
has been carried out using atmospheric neutrino data observed 
in Super-Kamiokande-I. The tau neutrino excess events have 
been observed in the upward-going direction as expected. The 
Super-Kamiokande-I atmospheric neutrino data are consistent 
with <-> v T oscillations. 
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Mining and Smelting Company. The Super-Kamiokande ex- 
periment has been built and operated from funding by the 
Japanese Ministry of Education, Culture, Sports, Science and 
Technology, the United States Department of Energy, and the 
U.S. National Science Foundation. 
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